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^ \ Abstract 

We constrain an interacting, holographic dark energy model, first proposed by two of us in |l|, with 
observational data from supernovae, CMB shift, baryon acoustic oscillations, x-rays, and the Hubble rate. 
The growth function for this model is also studied. The model fits the data reasonably well but still 
the conventional ACDM model fares better. Nevertheless, the holographic model greatly alleviates the 
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coincidence problem and shows compatibility at la confidence level with the age of the old quasar APM 
O ' 08279+5255. 
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I. INTRODUCTION 

Models of holographic dark energy have become popular as they rest on the very reasonable 
assumption that the entropy of every bounded region of the Universe, of size L, should not exceed 
the entropy of a Schwarzschild black hole of the same size. That is to say, 

L 3 A 3 < S B h =* L 2 M 2 Pl (M Pl = (SitG)~ l ) , (1) 

where A stands for the ultraviolet cutoff, the infrared cutoff is set by L. 

However, as demonstrated by Cohen et al. [2J, an effective field theory that saturates the above 
inequality necessarily includes states for which the Schwarzschild radius exceeds L. It is therefore 
natural to replace the said bound by another one that excludes such states right away, namely, 

L 3 A 4 < M Pl L . (2) 

This bound guarantees that the energy L 3 A 4 in a region of the size L does not exceed the energy 



of a black hole of the same size |3j . By saturating the inequality ([2]) and identifying A 4 with the 
density of holographic dark energy, px , it follows that [2j, |3 1 

3 C 2 

Px = mgT^ (3) 

where the factor 3 was introduced for convenience and c 2 is a dimensionless quantity, usually 
assumed constant, that collects the uncertainties of the theory (such as the number particle species 
and so on) . For a more thorough motivation of holographic dark energy see Section 3 of [l| . 

Last relationship is widely used in setting models of holographic dark energy that aim to explain 



the present stage of cosmic accelerated expansion, 



% 



7| , via the huge negative pressure associated to 



them. Broadly speaking holographic dark energy models fall into three main groups depending on 
the choice of the infrared cutoff, L. Namely, the Hubble radius Lll-. LSfl , the event horizon radius [3,^- 

n nn n 

Il3l ]. and the Ricci's length [l44l7l|. The particle horizon radius was also used |18[] but it presents 
the severe drawback of leading to a cosmology incompatible with a transition from deceleration to 
acceleration during the Universe expansion. 

In this paper we consider a spatially flat Friedmann- Robertson- Walker universe dominated by 
holographic dark energy (with the infrared cutoff set by the Hubble radius, i.e., L = H' 1 ) and 
pressureless dark matter such that these two components are dynamically linked by an interaction 
term. The model was introduced in [if]. Here we constrain it with data from supernovae type 
la (SN la), the shift of the first acoustic peak in the cosmic background radiation (CMB shift), 



baryon acoustic oscillations (BAO), and x-rays (strongly related to the baryon gas abundance in 
galaxy clusters), and Hubble's history, H(z). We also study the evolution of the growth function 
which potentially may constrain the model as well. But we do not use them because, at present, 
these data are far noisier than those in the other data sets. 

The paper is organized as follows. Section II recalls the holographic interacting model. Section 
III constrains the model with observational data. Notwithstanding it does not contains the flat 
ACDM model as a limiting case (at variance with, e.g., quintessence models) it shows a sizable 
overlap with the latter. Section IV studies the growth function. Finally, Section V summarizes 
our overall conclusions. As usual, a zero subindex indicates the present value of the corresponding 
quantity. 

II. BASICS OF THE MODEL 

n 

The spatially flat FRW holographic model proposed in [1[ rests on two main assumptions: (i) 
The dark energy density is governed by the saturated holographic relationship, Eq. ([3]), with the 
infrared cutoff fixed by the Hubble radius, i.e., L = H . (ii) Dark matter and dark energy do not 
evolve independently of each other. They interact according to 

p M + 3iTpM = Q , and p x + 3iT(l + w)p x = -Q , (4) 

where w = px/px stands for the equation of state parameter of dark energy, and 

Q = Tp x (5) 

is the interaction term where T denotes the rate by which px changes as a result of the interaction. 
We assume T to be semipositive-definite. Note that if Q were negative, the transfer of energy would 
go from dark matter to dark energy, in contradiction with the second law of thermodynamics [19]. 
Further, use of the Layzer-Irvine equation on nearly one hundred galaxy clusters strongly supports 
this view 201 ] . 



Interacting models were first proposed by Wetterich to lower down the value of the cosmological 



term 21]. Later on it was proved efficient in easing the cosmic coincidence problem 22j, [23| and it 



was suggested that the interaction (whatever form it might take) is not only likely but inevitable 



24, 



251 ] . The amount of literature on the subject is steadily increasing -see, e.g., [26j] and references 
therein. Admittedly, the expression (JSJ) is nothing but a useful parametrization of the interaction. 
Given our poor understanding of the nature of dark matter and dark energy, there is no clear 



guidance to derive an expression for Q from first principles. This is why our approach will be 
purely phenomenological. 

The model is fully specified by three quantities, e.g., the current value of the Hubble rate, Hq, 
the dimensionless density parameter Qx '■= & , xGpx/{3H 2 ) (or, equivalently, Qm)i and T. Note 
that c 2 is fixed by c 2 = Qx, as it can be readily checked. 

The first assumption readily implies that £lx does not vary with expansion, and that the ratio 
of energy densities, r := pm/px, stays fixed in spatially flat FRW universes (Qm + &>x = 1) for 
any interaction. The latter consequence greatly alleviates the coincidence problem albeit, strictly 
speaking, it does not solve it in full because the model cannot predict that r ~ 0(1) (to the best of 
our knowledge, no model is able to predict that). This feature of £lx and r being strictly constants 
may seem too strong; however, one should bear in mind that both quantities would slightly vary 
with the Universe expansion if the parameter c 2 in Eq. (J3j) were allow to weakly depend on time, 
something not at all unreasonable. Further, r would not be constant if the restriction to spatial 
flatness were relaxed. At any rate, we shall take the conservative stance that both c 2 and r do not 
vary; thus, the number of free parameters of the model will be kept to a minimum. 

At first sight, the consequence of Ox being of order unity also at early times might look wor- 
risome. One may think that a large dark energy component at that period would prevent the 
formation of gravitationally bound objects. However, this is not the case as 

i + r r 

W = -—3H (6) 

is not constant, and for suitable choice of the ratio T/H it tends to the equation of state of 
non-relativistic matter at early times. Its evolution is governed by the Hubble rate which, in the 
simplest case of V being a constant, takes the form 

»=*[*h + ( i -*ky-*']' < 7 » 

which corresponds to a specific generalized Chaplygin gas [271 ] . In last expression, the scale factor 
has been normalized by setting ao = 1. Figures Q] and [2] show the history of the equation of state 
for the best fit values of the model up to redshifts 8 and 1.2, respectively. Figure [T] illustrates that 
at high redshifts w approaches zero asymptotically. Figure [2] shows that, in accordance with the 



analysis in 28J, w(z) varies little at small redshifts. 

The deceleration parameter, q := —a/(aH 2 ), whose evolution is illustrated in Fig. [3l obeys 




Figure 1: Evolution of the equation of state parameter of dark energy, Eq. ([6]), for the best fit model, up to 
z = 8. In this, as well as in subsequent figures, the red swath indicates the region obtained by including the 
1<7 uncertainties of the constrained parameters used in the calculation (in the present case, VLx and T / Ho). 
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Figure 2: Evolution of the equation of state of dark energy for the best fit model up to z = 1.2. The 



observational data with their 2<r error bars are borrowed from 
was made. 



28] . In plotting the curve no fit to these data 



This expression implies that q — > * a t high redshifts as it should, and that the transition from 



deceleration to acceleration occurs at 



Ztr 



2r 



3# n r - T 



2/3 



1 



(9) 



which yields ztr — 0.80 for the best fit values. It should be noted that in [8| the transition 
deceleration-acceleration required that the c 2 varied, if only very slowly. In the present case, the 
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Figure 3: History of the deceleration parameter, according to Eq. |[8} in terms of redshift for the best 
fit holographic model (solid line). The redshift at which the transition deceleration-acceleration occurs is 
approximately 0.80. Also shown is the prediction of the ACDM model (dashed line). In this, as well as in 
subsequent figures, the green swath indicates the region obtained by including the la uncertainties of the 
constrained parameters used in the calculation (in the present case just fij^o)- The data are borrowed from 
|29j. In drawing the curves no fit to these data was made. 



transition also occurs for c 2 = constant (as, for simplicity, we are considering). The difference 
stems from the fact that in [8] the ratio T/H was kept constant, while the present model has T = 
constant, instead. 

The age of old luminous objects at high redshifts can constrain cosmological models by simply 
requiring that their age at the redshift they are observed do not exceed the age of the Universe 
at that redshift. Figure H] depicts the dependence of the age of the Universe on redshift for the 
best fit values of both the holographic model and the ACDM model alongside the age and redshift 
of three luminous old objects, namely: galaxies LBDS 53W069 (z = 1.43, t = 4.0 Gyr) [30] and 



LBDS 53W091 (z = 1.55, t = 3.5 Gyr) 



31 



32j, as well as the quasar APM 08279+5255 (z = 3.91, 



t = 2.1 Gyr) 33j, |34|. While the ages of the two first objects are lower than the ages of the 



holographic model and the ACDM model at the corresponding redshifts, the age of the quasar 
APM 08279+5255 lies slightly further than la beyond the age of the ACDM model at z = 3.91. 
By contrast, the holographic model is compatible at la level with the age of the said quasar. The 
tension between the APM quasar and the ACDM model has been known for some time now (see 



34| and references therein) and it has been revisited recently 351.136 1. 
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Figure 4: Dependence of the age of the Universe on redshift for the holographic model (solid line) and the 
ACDM model (dashed line). Also shown are the ages and redshifts of three old luminous objects, namely: 
galaxies LBDS 53W069, and LBDS 53091, and the quasar APM 08279+5255 - the latter with its ler error 
bar. In plotting the curves we have used the best fit value Hq = 68.1 ±2.1 km/s/Mpc for the holographic 
model and H = 72.1^ km/s/Mpc for the ACDM model. 



III. OBSERVATIONAL CONSTRAINTS 



In this section we constrain the three free parameters (fix, T/Hq, and Ho) of the holo- 
graphic model presented above with observational data from SN la (557 data points), the 
CMB-shift, BAO, and gas mass fractions in galaxy clusters as inferred from x-ray data (42 
data points), and the Hubble rate (15 data points) to obtain the best fit values. As the 
likelihood function is defined by C ex exp(— x 2 /2) the best fit follows from minimizing the sum 



Atotal 



Xsn ' Xcmb ' Xbao ' Xx—rays ' Xnubble' 



A. SN la 



We contrast the theoretical distance modulus 



Hth(zi) = 5 log 



10 



Dl 

lOpc 



+ Mo, 



(10) 



where ^o = 42.38 — 51og 10 /i, with the observed distance modulus fiobs(zi) of the 557 supernovae 
type la assembled in the Union2 compilation [7]. The latter data set is substantially richer than 
previous SN la compilations and presents other advantages; mainly, the refitting of all light curves 
with the SALT2 fitter and an upgraded control of systematic errors. In (fTOj) Dl = (l + #) Jq e(z'-p) 
is the Hubble-free luminosity distance, with p the model parameters {fix, T/Hq, and Ho), and 
E(z;p):=H(z;p)/H . 

The x 2 from the 557 SN la is given by 

2 I \ V^ i^thizi) - (J-obs\ z i)\ , . 

x sn (p) = ^ ^ , (11) 

where <j% stands for the la uncertainty associated to the ith data point. 

To eliminate the effect of the nuisance parameter /j,q, which is independent of the data points 
and the data set, we follow the procedure of [37] and obtain xl n = ^ s ^ nintmurn ' = 569.497. 

B. CMB shift 

The CMB shift parameter measures the displacement of the first acoustic peak of the CMB 
temperature spectrum with respect to the position it would occupy if the Universe were accurately 
described by the Einstein-de Sitter model. It is approximately model-independent and given by 



38, 



39] 



fZrec J- 

where z rec ~ 1089 is the redshift at the recombination epoch. The 7-year WMAP data yields 
K{zrec) = 1-725 ± 0.018 [a]. The best fit value of the model is K(z rec ) = 1.753±oQp. Minimization 

of 

XcmbiP) = 2 ( 13 ) 

a n 

-, 2 (minimum) o onr 

produces X CM B-shift = 2 - 385 - 



BAO 



Baryon acoustic oscillations can be traced to pressure waves at the recombination epoch gen- 
erated by cosmological perturbations in the primeval baryon-photon plasma. They have been 



revealed by a distinct peak in the large scale correlation function measured from the luminous red 
galaxies sample of the Sloan Digital Sky Survey (SDSS) at z = 0.35 [40J, as well as in the Two 
Degree Field Galaxy Redshift Survey (2dFGRS) at z = 0.2 [41[. The peaks can be associated 
to expanding spherical waves of baryonic perturbations. Each peak introduces a characteristic 
distance scale 



D v (zbao) 



ZBAO 



zbao dz 



H{z BAO ) \J H(z) 

(see Ref. [42j for a pedagogical derivation of this expression). 

Data from SDSS and 2dFGRS measurements yield A, (0.35)/ A, (0.2) = 1.736±0.065 



(14) 



41J] . The 



best fit value for the holographic model is A(0.35)/A (0.2) = 1.642 ± 0.003, and minimization of 



XbaoiP) 



2 (minimum) n „„„ 

g lves X b ao ; = 2.089. 



([A,(0.35)/A,(0.2)U - [A, (0.35) /A, (0.2)] 



obs ) 



(15) 



'D„(0.35)/D„(0.2) 



D. Gas mass fraction 



Since the bulk of baryons in galaxy clusters are in the form of hot x-ray emitting gas clouds (other 
baryon sources lagging far behind in mass) the fraction of baryons in clusters, f gas '■= M gas /Mtot, 
results of prime interest for it seems a good indicator of the overall cosmological ratio fibaryons/^M 



43J. The aforesaid fraction is related to the 



and, up to a fair extent, it is independent of redshift 

cosmological parameters through f gas oc d^ , where 

d,A '■= (1 + z)^ 1 Jq h(z') s * an ds for the angular diameter distance to the cluster. 

We used 42 Chandra measurements of dynamically relaxed galaxy clusters in the redshift interval 
0.05 < z < 0.1 



44]]. To fit the data we have employed the empirical formula 

JgasyZ) 



KA 7 b(z)n B0 ( djf DM \ 3/2 



(16) 



1 + s(z) n M o \ d A J 
(see Eq. (3) in Ref. [44]) in which the ACDM model is utilized as reference. Here, the parameters 
K, A, 7, b(z) and s(z) model the amount of gas in the clusters. We fix these parameters to their 
respective best fit values which can be found in Ref. [441 ] . 
The x 2 function from the 42 galaxy clusters reads 
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X 



x—rays 



(p) = E 



i=l 



\[Jgas\Zi)\th [Igas \%% ) jobs) 



(17) 



and its minimum value results to be Xx-""-™ 1 '" = 44.758. 



^-x-rays 



Figure [5] shows the fit to the data. 
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Figure 5: Gas mass fraction in 42 relaxed galaxy clusters vs. redshift. The solid and dashed curves 
correspond to the best fit models: holographic and ACDM, respectively. The data points with their error 



bars are taken from Table III in Ref. 



44] 



E. History of the Hubble parameter 



The history of the Hubble parameter, H(z), is poorly constrained though, recently, some 
high precision measurements by Riess et al. at z = 0, obtained from the observation of 



240 Cepheid variables of rather similar periods and metallicities 45], and Gaztahaga et al, at 
z = 0.24, 0.34, and 0.43 [46], who used the BAO peak position as a standard ruler in the radial 
direction, have improved matters somewhat. To constrain the model we have employed these four 
data alongside 11 less precise data, in the redshift interval 0.1 ^ z < 1.8, from Simon et al. [47J 



and Stern et al. 48J, derived from the differential ages of passive-evolving galaxies and archival 
data. 



Minimization of 



15 



XHubbleVP) — / , 



i=l 



[H th (zj) - H obs {zj)f 
a 2 {zi) 



(18) 



. , -, .,-, 2 (minimum) 

provided us with x Hu bble 



11.897 and H = 68.1 ± 2.1km/s/Mpc as the best fit for the 
Hubble's constant. Figure O depicts the Hubble history according to the best fit holographic model 
alongside the best ACDM model. 
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Figure 6: Plot of H(z) for the best fit values of the holographic model (solid line) and the ACDM model 
(dashed line). The data points and error bars are borrowed from Refs. 451-147 1. 



Figures [7] and [8] summarize our analysis. The left panel of Fig. [7] depicts the 68.3% and 95.4% 
confidence contours for SN la (orange), CMB shift (brown), BAO (blue), x-ray (black), and H(z) 
(green), in the (Qx, T/i^o) plane. The joined constraints corresponding to Xtotai are snown as 
shaded contours. The right panel depicts the 68.3% and 95.4% confidence regions in the (£lxo, Hq) 
plane of the holographic model (shaded regions) and the ACDM model (blue contours). As it is 
apparent, the models present a non-small overlap at 2a level. 

Figure [8] depicts the normalized likelihoods, £ oc exp(— Xtotai/^)' °^ * ne three free parameters of 
the holographic model. 

Altogether, by constraining the holographic model presented in Section II with SN la, CMB- 
shif, BAO, x-rays, and H(z) data we obtain Q x = 0.745 ± 0.007, T/H = 0.563±g:gJ|, and H = 
68.1 ± 2.1km/s/Mpc as best fit parameters, with Xtotai = 630.627. This value lies well inside the 
la interval {Xtotai/dof ~ 1-03). It should be noted that the no interacting case is discarded at 
very high confidence level. This means no surprise at all since for T = the model reduces to 
the Einstein-de Sitter (Qm = 1, &x = 0) and accordingly, as Eq. ([8]) tells us, the transition from 
deceleration to acceleration cannot occur. 

Table U shows the partial, total, and total x 2 over the number of degrees of freedom of the 
holographic model along with the corresponding values for the ACDM model. In the latter one 
has just two free parameters, SI mo an d Hq. Their best fit values after constraining the model to 
the data are D, M o = 0.259±g;gg|, and H = 72.1±J;| km/s/Mpc, with x\ Q tai = 593.142. 
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Figure 7: Left panel: the 68.3% and 95.4% confidence contours for the pair of free parameters ($lx , T/Ho) 
obtained by constraining the holographic model with SN Ia+CMB-shift+ BAO+x-ray+H(z) data. The 
joined constraints corresponding to Xtotal are rendered as shaded contours. The no interacting case is 
largely disfavored by the data. Right panel: the 68.3% and 95.4% confidence contours for the pair (fixo , 
Ho) of the holographic model (shaded contours) and the ACDM model (blue contours). The solid points 
signal the location of the best fit values. Notice the overlap at 2a confidence level between both models. 
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Figure 8: The normalized likelihoods of ilx, T/Hq, and H . 

We see that the ACDM model fits the data better than the holographic model in spite of 
having one parameter less. Thus, the former model should be preferred on statistical grounds. 
Nevertheless, this does not tell the whole story; the ACDM cannot address the cosmic coincidence 
problem and has some tension with the age of the old quasar APM 08279+5255. By contrast, the 
holographic model answers the said problem and shows compatibility, at la, with the age of the 
old quasar. 
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Model 


Xsn 


Acm6 


&bao 


X,x—rays 


X% 


Atotal 


XtotJdof 


Holographic 


569.497 


2.385 


2.089 


44.758 


11.897 


630.627 


1.03 


ACDM 


541.833 


0.013 


1.047 


41.527 


8.727 


593.142 


0.97 



Table I: x 2 values for the best fit holographic model (Cl x = 0.745 ± 0.007, T/H = 0.563±g$£, and 
H = 68.1 ± 2.1km/s/Mpc), and the best fit ACDM model (Q M o = 0.259^;^, and H = 72.1±^| 
km/s/Mpc). 

IV. EVOLUTION OF THE GROWTH FUNCTION 

It is not unfrequent to find in the literature cosmological models that differ greatly on their basic 
assumptions but, nevertheless, present a rather similar dynamical behavior. It is, therefore, rather 
hard to discriminate them at the background level. However, their differences are more readily 
manifested at the perturbative level (though, admittedly, the uncertainty in the corresponding 
data are, in general, wider). This justifies our interest in studying the evolution of the matter 
perturbations of the holographic model inside the horizon. 

A prime tool in this connection is the growth function, defined as 

f :=dlnS M /d\na, (19) 

where 8m denotes the density contrast of matter. In order to derive an evolution equation for /, 
we start from the energy balance for the matter component in the Newtonian approximation 

Om n v M = hi H • (20) 

a PM PM 

Here, vm is the velocity potential, defined by UMa = VM,a, where UMa is the matter four- velocity, 
and the hat means perturbation of the corresponding quantity. 

Recalling Eqs. (j3|) and © and that T and r do not vary, we can write 

k 2 r 

8m ?vm = — {hi - Sx) ■ (21) 

a z r 

Usually, the density contrast of dark energy is neglected under the assumption that dark energy 
does not cluster on small scales. However, as forcefully argued by Park et al. 49J], the neglecting 



of 5x can be fully justified in the case of the cosmological constant only. At any rate, in the 
present case the setting of 5x to zero wold be incorrect given the coupling between both energy 
components at the background level (i.e., Eqs. (j3])). It seems therefore reasonable to include 
a coupling, at least approximately, also at the perturbative level. The simplest possibility is to 

13 



assume a proportionality 5x = ol 5m with a constant a. As we shall see, the only consistent choice 
for this constant (under the conditions that T and r are held fixed) is a = 1. Thus, Eq. (|21|) 
becomes 



A; 2 T 
<5m n«M = (1 -a) 6m- 



(22) 



or r 

An equation for vm follows from the momentum conservation of the matter component. Assuming 
that there is no source term in the matter rest frame, this equation takes the simple form 



VM + <P = ; 



(23) 



where <f> is the Newtonian potential. Differentiation of (|22p. use of f|23|) and (|22p . and substitution 
of the scale factor for the time as independent variable, leads to 



X" i 



r 2(1 -«) r 



6\ 



M 



3 r + a 
2^2 r + 1 



1 



4(1 -a) r r + 1 
3 flV r + a 



5m = : 



(24) 



3Hr 3 flV 

where use of Friedmann's equation, 4irGp 
with respect to a. 

For a vanishing T we must recover the conventional perturbation equation 6 M -\-i^6' M — ^?5m = 
with the growing solution 6m °^ a f° r a dust universe. Clearly, this is only feasible for a = 1. With 
this choice the fractional matter perturbation 6m coincides with the total fractional energy density 
perturbation, 5 
holographic models reduces to 



pm px jt follows that the basic matter perturbation equation for the interacting 



6'L + 



'M 



2a 



1 + 



3Hr 



6', 



M 



2a- 



;6m = . 



(25) 



Replacing 6m by the growth function /, last equation becomes 



' ,+ ' 2+ K i+ iV-H 



(26) 



with /' := df/dlna. This has the advantage of being a first order differential equation. Notice that 
in the absence of interaction, T = 0, its solution is simply / = 1 as it should, i.e., a dust dominated 
universe. 

Figured] depicts the evolution of the growth function in terms of the redshift for the holographic 
as well as for the ACDM model. The latter appears to fit the data below z ~ 0.6 better than the 
former. In particular, at z = 0.15 the best fit holographic model deviates A/ = 0.3 (corresponding 
to 3cr) from the observed value (though it falls within la with the remaining data points) while 
the best fit ACDM model falls within la also at z = 0.15. 
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At any rate, it has been recently pointed out, from the observation of nearby gala xies, that 



5lf]. Clearly 



structure formation must have proceed faster than predicted by the ACDM model 

slightly enhanced values of / at low redshifts helps accelerate the formation of galaxies and clusters 

thereof. 
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Figure 9: Growth function vs. redshift for the best fit holographic model (solid line). Also shown is the 
jrediction of the ACDM model (dashed line). The observational data are borrowed from Table II in Ref. 
50| . In plotting the curves no fit to these data was made. 



V. CONCLUDING REMARKS 

We constrained the interacting holographic model of Section II with data from SN la, CMB 
shift, BAO, the gas mass fraction in galaxy clusters, and H(z). The parameters of the best fit 
model are: n x = 0.745 ±0.007, T/H = 0.563±£^ , and H = 68.1 ± 2.1 km/s/Mpc. We have not 
included data of the growth function in the likelihood analysis given the wide uncertainties of the 
current data. However, we have derived the differential equation for /, Eq. (I26p . and integrated it 
numerically for the best fit model. 

It conforms reasonably well to the observational data but not so well as the ACDM model 
(best fit values: Q M o = 0.259+o;°o5, H = 72.lt\i km/s/Mpc) does notwithstanding the latter 
has one less free parameter than the former. However, the holographic model greatly alleviates 
the cosmic coincidence problem and seems compatible at la level with the age of the old quasar 
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APM 08279+5255. Besides, the observational data from the CMB shift, BAO, x-ray, and some 
of the H(z) data, are not fully model independent owing to the fact that they are extracted with 
the help of the conventional ACDM. This frequently makes the latter tend to be observationally 
favored over any other cosmological model. Moreover, the BAO data are conventionally determined 
under the assumption of purely adiabatic perturbations. However, as recently argued [52|, should 
isocurvature components be present the shape and location of the CMB acoustic peaks would be 
altered and the data extracted from BAO affected. 

Clearly, we must wait for more abundant, varied, and model- independent accurate data to tell 
which of the two models survives. If eventually neither of the two does, we should not be so much 
disenchanted because, at any rate, this "negative" result would have narrowed significantly the 
parameter space of dark energy. 
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